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.EVALUATION OF UPWELLING INFRARED 
RADIANCE FROM EARTH'S TROPOSPHERE 


by 

S, K. Gupta and S. N, Tiwari 
School of Engineering 
Old Dominion University 
Norfolk, Virginia 23508 


SUMMARY 

Basic equations for calculating the upwelling atmospheric radiation 
are presented which account for various sources of radiation coining out at 
the top of the atmosphere. The theoretical formulation of the transmittance 
models (line-by-line and quasi-random band model) and the computational ’ 
procedures used for the evaluation of the transmittance and radiance are 
discussed in detail. By employing the Lor'entz line-by-line and quasi-random 
computer programs, model calculations were made to determine the upwelling 
radiance and signal change in the wave number interval of CO fundamental 
band. These results are useful in determining the effects of different 
interfering molecules, water vapor profiles, ground tanperatures , and ground 
emittances on the upwelling radiance and signal change. This information is 
of vital importance in establishing the feasibility of measuring the concen- 
trations of pollutants in the atmosphere from a gas filter correlation instru- 


ment flown on an aircraft or mounted on a satellite 



. • I, INTRODUCTION' 


The present status of knowledge about the fate of the relevant 
pollutants in the atmosphere Is limited. Determination of the concentra- 
tion of various atmospheric pollutants is important because of their con- 
siderable influence on the weather and environment. 

The earth's surface with its temperature in the vicinity of 300®K 
emits like a black body from the near to the far Infrared region of the 
spectrum. The emission in the near infrared (between 2 and 20 microns) is 
particularly important because most of the minor atmospheric constituents 
(i.e, , CO 2 , N 2 OJ H 2 O, CO, CH^, NH 3 , etc.) absorb and emit in this spectral 
region. The upwelling infrared radiation from the earth's atmosphere, 
therefore, consists of the modulated surface radiation and the radiation 
from the atmosphere. This radiation carries the spectral signature of all 
the minor atmospheric constituents amongst which gases such as CO, CH^ and 
NH 3 are called the atmospheric pollutants. 

Detection of the atmospheric pollutants from infrared radiation 
measurements' requires a detailed knowledge of the atmospheric constituents 
which absorb and emit significantly in the spectral range of Interest. A 
systematic representation of the absorption and emission by a gas, in the 
infrared, requires the identification of the major infrared bands and eval- 
uation of the line parameters of these bands. The line parameters depend 
upon the temperature, pressure, and concentration of the absorbing molecules 
and, in general, these quantities vary continuously along a nonhomogeneous 
path through the atmosphere. With the availability of high resolution 
spectrometers, it is now possible to determine the line positions, intensi- 
ties, and half-widths of spectral lines quite accurately. As a result, the 



line parameters of the strong infrared bands of most minor atmospheric 
constituents are now known quite well [1-3], 

An accurate model for the spectral absorption coefficient is of 
vital importance in the calculations of the transmittance of an atmospheric 
constituent and in the correct formulation of the radiative flux equations 
which are employed in the reduction of data obtained from either direct or 
remote measurements. High spectral resolution measurements make it neces- 
sary to employ line-by-line models for data reduction. If, however, the 
integrated signals over a relatively wide spectral interval are measured, 
then one could employ an appropriate band model. The line models usually 
employed in the study of atmospheric radiation are Lorentz, Doppler, and 
combined Lorentz-Doppler (Voigt) line profiles. A complete formulation 
(and comparison) of the transmittance (and absorptance) by these lines, in 
an infinite and finite spectral interval, is given in [4-6]. The most 
appropriate band model for atmospheric application is the quasi-random 
model which is discussed in detail in this report, 

Ludwig et al. [3] have explored the possibilities of measuring the 
amount of atmospheric pollutants through- remote sensing. An important 
method of measuring the pollutant concentration by remote sensing is the 
passive mode (also called the nadir experiment) in which the earth-oriented 
detector receives the upwelling atmospheric radiation. The near infrared 
region is particularly suitable for passive mode measurements simply because 
the radiation in this region is practically free from the scattering effects. 
Radiation in the visible and ultra-violet regions is severely affected by 
the scattering processes which make meaningful passive mode measurements 
impossible. 



The purpose of .this study is to present a comprehensive analysis 
procedure for calculating the atmospheric transmittance and upwelling 
radiance in the wave nximber interval of CO fundamental band (between 
2070 and 2220 For this purpose computer programs were developed 

by employing the .Lorentz line-by-line* model as well as the quasi-random, 
narrow band model • The basic equations for calculating the upwelling 
atmospheric radiance are presented in Sec* 11* The theoretical formula- 
tions of the transmittance by the two spectral models is given in Sec* III. 
Listings of the computer programs are provided in the Appendicies* By 
employing the line-by-line and quasi-random computer programs ^ model cal- 
culations were made to study the effects of different interfering molecules, 
water vapor profiles, ground temperatures, 'and ground emittances on the 
upwelling radiance* These are presented in Sec* IV* 
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• II.- -UPWELLING- ATMOSPHERIC RADIATION •’ 

As shown in Fig, 2.1, the radiation. emergent from the atmosphere, 

E((i)) , may be given by the expression [3,5] 

E(w) = Eg(tc) + Ej^(w) + E^(U) + (2.1) 

where 

■ Eg(w) = thermal radiation emitted by underlying surface and 
atmosphere 

Eg^(to) = incident solar radiation reflected, by the surface 
E^((o) = radiation scattered by single or multiple scattering 

processes in the atmosphere without having been reflected 
from the surface 

E^^(w) = scattered energy which has undergone a reflection from 

the surface. 

In general, .these quantities are functions of surface temperature, atmo- 
spheric temperature, surface emittance, surface reflectance, , sun zenith 
angle, scattering characteristics of particles, and transmittance of the- 
atmosphere. 

In the spectral region of infrared measurements, the effect of scat- 
tering is negligible. The incident solar radiation reflected by the 
surface, however, is important especially if the surface reflectance is 
assumed to be high (0.2 and higher). 

Upon neglecting the scattering and solar radiation, the expression for 
thermal radiation emerging from a plane-parallel atmosphere can be written 
as 



Pollutants Atmosphere 



Fig. 2.1 Radiative energy received by an aircraft or satellite mounted instrument 
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E(o)) = E (to) - e(to) B.(to,T ) T(to,o) 
v s 


f 


B(to,T(z)) [dx(to,z)/dzldz 


( 2 . 2 ) 


where e(o)) is the surface emittance, B(to,T) is the Plank's blackbody 
I 

function, T is the surface temperature, T(z) is the temperature at 

s 

altitude z ,.and x(( 0 ,z) is the monochromatic transmittance of the 
atmosphere. The first term on the right hand side of this equation 
represents the radiation from the surf ace. while the second term is the 
radiation from the atmosphere. 

The contribution from sunlight reflected from the surface becomes 

significant at shorter wavelengths. This contribution is given by the 

component E^(o3) as 
K 

Ej^(to) = (1 /it)[ 1 - e((iJ)] cos e Hg(w) [x(o))]^ (2-.3) 


wher,e 6 is the sun zenith angle and ^ - 1 + f(0) * Function f (6) 

sec 0 for o ^ 0 ^ 60^ and equals to Ch 0 for 0 > 60®, with Ch 0 

denoting the Chapman function. H (o)) is the sun irradiance on top of 

s 

the atmosphere-, and x(to) ■ is the transmission vertically through the 
atmosphere. 

The expression for the transmittance x(w,z) is obtained by employ- 
ing an appropriate molecular (band or line) model, and, in general, this 
can be written as 


X(w,z) 


n [x (w,z)3 
i ^ 


exp[- 


exp[- 



2 e.(o),z') n.(z') dz*] 

• IL. 


2 3^(io,z') C^(z') P^(z’) dz»] 


(2.4) 
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where ’3^(uj,z) is the extinction coefficient for species i (which is 
equal to the sum of absorption and scattering coefficients), n^(z) is 
the number density of the contributing molecules of species i in the 
optical path, ^^( 2 ) is the concentration of species i , and P^(z) is 
the total pressure. If scattering is neglected, then atmospheric trans- 
mittance can be expressed in a convenient form as 


T(to,z) = niT^(aj,z)] = exp[- f '2 K“(aj,x*) p.(x’) dx^] 


(2.5) 


where ic™(oJ,x’) is the mass absorption coefficient for the ith absorbing 
gas, is the mass density, and x is the depth of the level measured 

from the top of the atmosphere. 

Upon introducing new definitions for the absorption coefficient and 
path length as [7,8] 

k(co) = fC™((o) (Pj./^g) 9 ( 2 * 6 ) 



(2.7) 


( 2 . 8 ) 


where P is the partial pressure and p is the mass density of the 
8 S 

absorber, u is the optical path length in cm atm, and p^ is the 
absorber density corresponding to reference conditions (usually chosen 
as STP) . 

A combination of Eqs. (2.2) and (2,8) yields a relation for thermal radia- 
tion emerging from a plane-parallel atmosphere. This equation, with appro- 
priate spectroscopic information, is used to obtain the concentration of 



atmospheric pollutants from, radiation measurement’s provided other governing 
parameters are knovm. 

In actual calculations of atmospheric transmittance, the nonhomogeneous 
atmosphere is' divided into a number of layers and for each layer an average 
value of pressure and temperature is assumed. With reference to Eq, (2.7), 
the pressure path length for, each layer can be written as 




(2.9) 


where Q. . is the volume , mixing ratio of the ith constituent in the jth 
xj 

layer,' dx^ xs the thickness of the jth layer, and and are the 

average pressure and temperature of the jth layer* ■ Since the line Intensity 
depends upon temperature and the line-width on temperature and pressure both, 

if 

the absorption coefficient will be slightly different for different. layers. 
Equation (2.8) can, therefore, be written as 

►u 


T(to, z) = exp E- 


L ' ' 


(dj) du..] . 


( 2 -. 10 ) 


where ^ indicates the summation over all the layers between the top and 
level z . . . . 

In general, the absorption coefficient for a single line centered at 
the wave number is given by 


= S f (w, Y ) , 
con n n n 


( 2 . 11 ) 


where is the intensity of the nth spectral line and is given by 


S = /. K d(co-ra) ) 

n I con n . 

rvi ^ 


( 2 . 12 ) 
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The line intensity may be described in terms of the molecular number 
density and Einstein coefficients, i.e,, it depends upon the transition 
probabilities between the initial and final states and upon the popu- 
lations of these states. For a perfect gas it may be shown that 
is a function solely of temperature. The quantity Yjj) 

line shape factor for the nth spectral line. It is a function of the wave 

number (j) and the line half width Y and is normalized on (w-O) ) 

n ' n 


such that 


f f (w-o) ■ 

I n n 

CO 


) d(w-u^) = 1 


Several approximate line profiles have been described in the 
literature. Most commonly used profiles are rectangular, triangular, 
Lorentz, Doppler, or Voigt (combined Lorentz and Doppler) profiles. 
Lorentz, Doppler, and Voigt profiles are of special interest in the 
atmospheric studies and these are discussed in detail in [4,6]. 

As indicated in Refs. [4,5,9], use of the Lorentz line profile is 
justified for radiative transfer analyses in earth's lower atmosphere. The. 
shape factor, f^ , for the nth rotational Lorentz line is such that the 
absorption coefficient for species i can be written from Eq. (2.11) as 


(k. .) = S f (w,Y ) = S- Y /{Tr[(o>-w )^ + Y^]} . 

ij n n n * 'n n 'n ^ n'^ 'n 


It should be emphasized here that the units for the absorption 
coefficient and the pressure path length, in Eqs, (2,6) and (2.7), are 
cm ^ - atm ^ and cm-atm respectively. In the case of water vapor, 
however, the path length is expressed in precipitable centimeters (pr - cm) a 
and the corresponding absorption coefficient in (pr - cm) 
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The total atmospheric transmittance is obtained by combining Bqs,(2*10) 
and (2*11) and integrating the resulting expression over the spectral range 
of interest. The exact procedure for calculating this, by employing the 
line-by-line and quasi-random band models, is discussed in the next section. 

The total energy emergent from the atmosphere is obtained by integrating 
either Eq,(2,l) or (2,2) over the specified spectral interval A(0 as 

j[^ E(m) du (2.15) 

The procedure for calculating the upwelling radiance, by employing the line- 
by-line and quasi-random band model for the transmittance, is briefly discussed 
in Sec. IV. 



III. EVALUATION OF ATMOSPHERIC TRANSMITTANCE 


As pointed out earlier, an accurate model for the spectral absorption 
coefficient is essential in the -calculation of atmospheric transmittance and 
in the correct formulation of the radiative flux equations which are used in 
the reduction of data obtained from either direct or remote measurements. 

While use of line-by-line models is essential in the data reduction of high 
spectral resolution measurements, an appropriate band model could be employed 
if integrated signals are measured over relatively wide spectral intervals. 
Several line-by-line and band models are available in the literature. In 
-this section, the procedure for obtaining the atmospheric transmittance by 
employing the line-by-line model and the quasi-random band model is discussed, 

3.1 Direct Integration (Line-By-Line Model) 

A computer program has been developed to determine the atmospheric 
transmittance by employing the Lorentz line-by-line model. The procedure 
for computing the transmittance is presented here and a listing of the 
program is given in Appendix B. Explanation of the symbols used in the 
computer program is given in Appendix A. The procedure can easily be 
extended to incorporate other line profiles. 

The direct integration method consists of calculating the absorption 
coefficient, and then the transmittance, at a large number of frequencies 
within the spectral range of interest. Since the absorption coefficient is 
a highly varying function of the frequency (varying by orders of magnitude 
over the width of a single line, which is often less than 1 cm ^) , it has 
to be evaluated at very closely spaced locations. The total absorption 
coefficient at any frequency location is made up of contributions from a 
large number of lines in the vicinity of that frequency. Some molecules 



have thousands of lines in the spectral range of a particular hand and^ the 
contribution from each one of these lines should be evaluated at each 
frequency location. It is, therefore, evident that the evaluation of the 
atmospheric transmittance with high accuracy is a time consuming task# 

Since the absorption coefficient varies rapidly only in the vicinity of 
the line center, it should be evaluated at very close intervals near the 
line centers, while away from the line centers bigger intervals could be 
used without any serious loss of accuracy# Kunde and Maguire [8] and 
Drayson [9] have proposed a scheme for selecting the frequency locations 
for the calculation of absorption coefficients and transmittances# The 
scheme essentially consists of numerically evaluating the average trans^ 
mittance over a narrow spectral interval by employing the Lang endre— Gauss 
quadrature formula [10]. The scheme employed in the present study is a 
modified form of the formulation presented in [8,9] where a different set 
of parameters are used. The choice of the parameters and the basic features 
of the computational scheme is discussed in the following subsections# 

3.1#1 Spectral Divisions into Intervals and Sub intervals 

The entire frequency range of interest (usually the spectral Interval 
of a particular vibration-rotation band) is first divided into a large 
number of rather narrow intervals Ato # Each of these Ao) intervals is 
then divided into a variable number k of sub intervals depending upon the 
number of lines falling in that interval# The arrangements of spectral 
divisions are illustrated in Fig. 3.1# Two very narrow subintervals are 
created on each side of the line center and these are denoted in the computer 
program by ALX (see Appendix B) . As shown in Fig. 3.1, an interval Aco 
remains one subinterval if no line is present, but is divided into 6 or 11 
sub intervals in the presence of one or two lines respectively. For greater 
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(a) 

NO LINE 
( 4 mesh points ) 


(b) 

ONE LINE 
( 24 mesh points ) 



i 



(c) 

TWO LINES 
( 44 mesh points ) 


M 


a narrow interval Ato o 


number of lines in the interval Au , the number of subintervals will 
also be greater.- Four frequency locations are then obtained within each 
subinterval by using the relation 

“it ■ "V - 

where k refers to the kth subinterval and are the upper and 

XXIX 

lower frequency liioits of the sub interval respectively, and ^ are the 

X/ 

abscissa values for the four point Gauss-Legendre quadrature formula [10], 

3.1.2 Evaluation of Absorption Coefficients 

The resultant magnitude of the absorption coefficient at the wave number 
location is due to two separate contributions; (a) direct contribution 

and (b) wing contribution and is given by the equation 

The contribution to the absorption coefficient from lines whose centers 
lie in close vicinity (on both sides) of the wave number under consideration 
is called the direct contribution and is obtained for Lorentz lines from 
Eq. (2.14) as 

^ ■ I (3.3) 

where o)^ refers to the center of the nth contributing line. In the computer 
program, the range of direct contribution is denoted by DLIM, Thus, if 
I (00-03^) [ ^ DLIM , the contribution is called the direct contribution and is 
evaluated by using Eq. (3.3). 

The wing contribution to arises from lines located at wave numbers 

farther than DLIM from ai „ (on both sides). The wave number range of the 
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wing contribution is denoted by WLIM and the .absorption coefficient for 
this range is calculated from the equation 

“ 2 = S {S Y /[ir(to-w y2]> , (3,4) 

n n ^ 

only for the interval DLIM < | (oj-o)^) ] ^ WLIM . The contribution from the 
lines for which | | > WLIM is found to be negligible. The range of 

WLIM is established by several trial calculations. 

Since the number of lines between DLIM and WLIM (on each side) is 
very large, the wing contribution is ’not evaluated at each individually. 

Instead, it is evaluated at the boundaries' and the center of an interval. The 
values at all the within that interval is obtained by linear interpo- 

lation. As pointed out by Drayson [9], this procedure does not affect the 
accuracy of the results seriously. 

Since varies with temperature and varies with temperature as 

well as pressure, the nonhomogeneous atmosphere is divided into appropriate 
numbers of homogeneous slabs for the purpose of calculating the absorption' 
coefficients and the atmospheric transmittances . For each slab (or layer) • 
an average value of the temperature and pressure is assumed. 


3.1.3 Optical Thickness and Transmittance 

The path length of each homogeneous slab of a nonhomogeneous atmosphere 
is evaluated separately. For any constituent, the path length for the jth 
layer is found from either Eq. (2.7) or Eq. (2.9) to be 


u. 

3 






(3.5) 


where is the thickness of the jth layer* 


The optical thickness of the jth layer at the wave number is 


given by 



“j •• 


(3.6) 


and tixe total optical thickness for the atmosphere by 

2 • . (3-7) 

For several absorbing species in the spectral range of interest, the total 
optical thickness is given by 

2 2 [k(01j^)_ , (3,8) 

where the summation over i constitutes the contribution of, all the absorb- 
ing constituents. 

The atmospheric transmittance at can now be expressed by 

= «qp{-S s Uy-1> , (3.9) 

where ic(a) ) is given by Eq. (3.2), 

KX 

The average, transmittance over the interval Au, is obtained from the 
equation 

t(4m) ^ T(03^p]}/4M , (3.10) 

where the summation k extends over all the subintervals within the interval 
A(0 , and weight factors for the four point Gauss-Legendre 

quadrature formula [10] . 

3.1.4 Optimization of Spectral Parameters 

The spectral parameters (such as Alo , ALX , DLIM , WLIM , .etc.) 
entering in the numerical calculation of atmospheric transmittance are 
optimized for accuracy and for saving the computational time,. Various 



parameters are carefully chosen on the basis of considerable numerical 
experimentation. In a particular numerical scheme for calculation; of the 
atmospheric transmittance^ the choice of an optimized spectral parameter 
may slightly depend on the nature of the- molecules under consideration 
(i.e. on the spectral distribution of the rotational lines). For the 
present study, spectral parameters were optimized by considering the (X) 
spectral lines in the region 2070-2220 cm ^ and by assuming a uniform CO 
distribution of one part per million (1 ppm) by voltime in the atmosphere 
from zero to ten kilometers. 

As pointed -out earlier, the entire spectral range of interest is- first 
divided into a large number of narrow intervals Aw’s . In the computer 
program, this parameter is denoted by DEL (see Appendix B) . By employing 
■Eq. (3.10) , the transmittance values were calculated for Aw “ DEL = 

2.0, 1.0, and 0'i5 cm**^ . These are given in Table 3.1 . It, is seen that 
while the transmittance increases by about 0.1% when DEL decreases from 
■2.0 to 1.0 cm ^ , it decreases. by about 0.03% when DEL goes from 1.0 to 
0.5' cm ^ . This variation being significantly small, the intermediate 
value of AciJ = DEL = 1.0 cm”^ is chosen for the present sttidy as -opposed 
to Ihe value of Ao) = 0.1 cm ^ used by Kunde and Maguire [8] . 

The effect of varying the width of the subintervals near the line centers 
(denoted by ALX in the computer program) was also studied by- numerical 
experimentation. • Average transmittance values were calculated for 
ALX = 0.04, 0.02, and 0.01 cm ^ and these are listed in Table 3.1. Since 
there is no effect of changing the ALX values on the transmittance, a value 
of ALX = 0.03 cm ^ was selected simply because the subintervals on- each side 
of the line center will cover a relatively larger portion of the important 
part of the line profile. 
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The effect of changing the value of DlilM on the average transmittance 
was found to be insignificant (see Table 3.1) and a value of DLIM = 5.5 cm ^ 
was selected for the present study in order to, save computational time. 

The effect of changing the WLIM values was, however, found to be significant 
and a value of WLIM = 45.5 cm"*^ was adapted. This value of \flLIM maintains the 
general accuracy of the program and provides savings of computational time. 

It was not necessary to consider a TOIM value beyond 45.5 cm ^ in view of the 
expected small gain in accuracy* 

• TABLE 3.1 

Average Atmospheric Transmittance from CO 
Fundamental Band as a Function of Various 
Spectral Parameters. 


Spectral 

Spectral 

Average Transmittance 

Parameter 

Parameter Range 
(cm“^) 

(nondimens ional) 



2.0 

0.79605 

DEL 

1.0 

0.79678 


0.5 

0.79655 


0.04 

0.79678 

ALX 

0.02 

0.79678 


0.01 

0.79678 

DLIM 

5.5 

0.79678 


7.5 

0.79679 


25.5 

0.79678 

WLIM 

35.5 

0.79655 


45.5 

0.79644 
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In the literature [1-3], individual ’line' widths ,are available for CO 

I 

and H 2 O lines while only the average values are available for CO 2 , N 2 O, 
and 0^ lines. For CO and H^O, model calculations indicate that the mean 
transmittances obtained by using the average line widths are different than 
those obtained by using the individual line widths. For CO, a value of 
average transmittance of 0,83282 was obtained by using the average line 
widths while a value of 0.82834 was obtained when individual line widths 
were used. Similarly for H 2 O, a value of T of 0.85563 was obtained for 
average line widths while a value of 0.86948 was obtained for individual 
line widths. In the. results of model calculations (reported in the next 
section), individual line widths were used for CO and H 2 O lines. 

The computer program listed in Appendix B calculates the atmospheric 
transmittance in the spectral 2070-2220 cm This is the range of CO 
fundamental band. The program accounts for the contribution of the inter- 
fering molecules such as H 2 O, CO 2 , N 2 O, and 0^. The program can easily be 
modified to calculate the atmospheric transmittance in any spectral range. 



3.2 Quasi-Random Band Model 

The use of narrow band models (such -as Elsasser, statistical, random- 
Elsasser, and quasi-random) Is' justified if integrated signals are measured 
over relatively wide spectral intervals. The use of these models in 
•transmittance calculations results in a considerable reduction in. compute-, 
tional time. The results obtained from most of these models, however, 

.usually lack the high accuracy needed in atmospheric transmittance calcula- 
tions, The quasi— random model, introduced by Wyatt j Stull, and Plass [11], 
is probably the best model- to represent the absorption of' a vibration- 
rotation band quite accurately. Discussions on various band model's are 
available in the literature [12-15] . The procedure for calculating the 
atmospheric transmittance by employing the quasi-random band model is dis- 
cussed here in detail. A listing of the computer program is given in 
Appendix C and various symbols used -in the program are explained in Appendix A 

The fundamental features of the quasi-random band model are discussed in 
references [7,11], For this model, the entire band span, A , is divided into 
a number of small subintervals of equal spectral width 6 .. The average 
transmittance is first evaluated for each 5 subinterval. The lines- within 
each sub interval are assumed to be distributed randomly. The average trans- 
mittance over the entire spectral range (i. e. , over the band pass) is 
obtained- by arithmetically averaging of the transmittances of all subintervals 

The lines in a particular subinterval 6 are divided into intensity 
subgroups. Reckoning from the intensity of the strongest line, five intensity 
groups (each spanning one decade) are created. Thus, lines whose intensities 
are within 10”® times the intensity of the strongest line (in that' group) are 
taken into consideration. It was suggested in [11] that the contribution 
from still weaker lines is negligible. 
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Since the line positions and intensities are asstimed to vary within 
the narrow subintervals 6 (rather than the entire band pass A ) , the 
results obtained from this model are expected to be better than other band 
models . 

3.2.1 Evaluation of Average ‘Transmittance - 

For a homogeneous path, the spectral transmittance at (0 due to a. ' 
line with center at is obtained by combining Eqs, (2,'8) and (2.11) as 

t(co) = exp[-S^ u = exp[-S^ u f(u,Uj^)] (3 


where the shape -factor f(co,{»)^) for the Lorentz line profile is defined 
in Eq. (2.14), 

For theoretical line shapes (i.e., for Lorentz, Doppler, and Voigt 
lines) the location of (O is symmetrical with respect to . Thus, 'the 
average transmittance of a line over 6 can be evaluated in two ways; 

(a) keep at its fixed location and average T(ti)) for location of w* 

varying over 6 , (b) keep ai fixed at the center of the subinterval 
and average t(o3) for location of .varying over 6 . The second approach 

is adapted in the calculation of- the average transmittance while .using the 
quasi-random band model. 

The. average transmittance over 6 due to a single spectral line may, 
therefore, be expressed by 


T (to) 
n 



u f(to,(o )] dto 
* n n 


> 


(3.12) 


where it should be noted that to^ is the variable of integration. 

If K is the niuBber of lines in an intensity decade then the average 
transmittance due to all the lines in that decade, , is given by 



22 


Td(w) =-£-| exp[-S^ u f(w,w^)] doj^}^ , 


where is the average intensity of all the lines within the decade 

under consideration. It is assumed that the range of intensity in a decade 
is not large enough to cause any serious errors because of intensity 
averaging, A considerable saving in computational time is realized if 
the number of lines in each intensity decade is large. 

The average transmittance due to all lines in the five intensity decades 
of the subinterval 6 is given by 

It 

where subscript k represents the.kth spectral subinterval (l.e. , 5^^) 

of the total spectral interval A , 

Equation (3.13) represents the transmittance due to the lines within 
, The wings of the lines in the adjacent subintervals also make a 
significant contribution to the absorption in , The resultant trans- 
mittance over the subinterval 6^^ , therefore, is given by 


n Vj - 


(3.14) 


where T, . (oi) represents the transmittance in 6. due to lines in 6 , . 

K-J K. 2 

Large variations are observed between the transmittances of adjacent 
subintervals if the above procedure is employed directly. Wyatt, Stull 
and Plass [11] attributed this to the arbitrary locations of the sub interval 
boundaries and the resulting arbitrary distribution of lines between the 
sub intervals . This spurious variation in transmittances is smoothed out by 
adapting the following procedure, (a) Divide the total spectral range A 
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(starting at A^) into K subintervalV with boundaries at + 5 , 

1 3 

A + 26 , A k6 and centers at A +*^,A etc* Evaluate 

o o2^o2’ 

1 3 

transmittances at the centers* i*e*at A +-;^,A etc* 

.02 o Z 

(b) Shift the spectral mesh of subintervals- by= • fi/2- The ,subinterv.al 
centers now lie at A , A. + 5 , A + ,26 , etc. Evaluate, transmittances' ’ 

O O'*’ » Q, - , . ^ f . 

at these shifted spectral locations. The resultant transmittances at the 
centers of original set of subintervals are obtained from these two sets by 
using the three point- sliding average procedure. For subinterval between 
A and A + 6 , the resultant transmittance at the center A -A +6/2 
is given by 

The' average transmittance for the entire range A is now expressed by 


t(A) = A 




K 

2 


k=l 




K 

2 

k=l 




(3.16) 


From the above description, it is evident that the basic quantity to 
be evaluated is the average transmittance due to a single line over the 
subinterval 6 . This is given by Eq. (3.12), Upon introducing the follow- 
ing change of variables 


y = 

1 

s 

o 

- 6/2 , 

(a) 

z = 

0 ) - to • - 
0 

6/2 , 

(b) 

'n 

S u/ (tt 
n 


(c) 

P = 

2y„/6 . 


(d) 

ri = 

2y/6 , 


(e) 

e « 

2z/6 , 


(f) 


(3.17) 
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Eq. (3.12) can.be expressed for the Lorentz line shape as 

= ( 1 / 2 ) f exp{-p^ ?„/[(e-Ti)^ + .p^]> dn . 


n 


(3.18) 


For direct contribution (i.e., contribution to the transmittance due 
to lines in the. same subinterval), e = o and the integrand in Eq. (3.18) 
is symmetrical about '0 = 0 . For this case, Eq. (3.18) reduces to 


■I 


T^(oj) = I exp[-p2 £ /(n^ + p^)] dn , 


(3.19) 


For evaluating the wing contribution, it is assximed that p^«(e-n)^ » 
and Eq. (3.18) reduces -to 


- (1/2) expi-p2 ^ 

The numerical procedures used for evaluating the integrals in Eqs. 
(3.19) and (3.20) are discussed in the next subsection. 


(3.20) 


3.2.2 Optimization of Integration Procedure and Parameters 

The integrand in Eq. (3.19) is a strong function of 0 and, therefore, 
an accurate procedure for evaluating the integral is required. Wyatt, Stull' 
and Plass [11] have obtained an analytical solution of Eq. (3.19), Following 
the procedure sugg^ted by Young [16], Kunde [7] divided the range of 
integration (0-1) into seven subintervals as 


.ri = 0.0, 0.001, 0.005, 0.01, 0.05, 0.1, 0.5, 1.0 , 

and used an eight point Gauss-Legendre quadrature formula for each subinterval. 

In the present work, a much simpler scheme was used to evaluate Eq.(3.19). 
Twenty-six abscissa points were chosen in the interval (0-1) -and suitable 
subintervals around them were created by a graphical experiment. Table 3.2 
lists the abscissa values and the subintervals associated with them. The 
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integrand is evaluated* at the twenty-six abscissa points, multiplied by the 
corresponding sub interval width, and the results are then summed. A com- 
parison of results of this procedure with the results of ten subinterval 
ten-point Gauss -Legendre quadrature formula (for the same interval) , 
revealed excellent agreement. It is, therefore, concluded that this , pro- 
cedure yields results of high accuracy. 

The integrand in Eq. (3.20) is a slow varying function of r| . Thus, 
a simple twenty-one point Simps on rule quadrature formula was used to 
evaluate the integral instead* of a ten-point Gauss-Legendre quadrature 
formula suggested by Kunde [7]. Accuracy checks were made on this quadrature 
scheme, and no appreciable differences were observed. 

, It should be emphasized here that the. use of these simple quadrature 
schemes results in considerable reduction in computational time. 

TABLE 3.2 


Abscissa Points and Sub intervals Selected 
to Evaluate the Direct Contribution. 


Abscissa 

Subinterval 

Abscissa 

Subinterval 

0.0 

0.0006 

0.04 

0.01 

0.001 

0.0006 

0.05 

0.01 

0.0015 

0.0006 

0.06 

0.015 

0.002 

0.0007 

0.08' 

0.02 

0.003 

0.001 

0.10 

0.03 

0.004 

0.001 

0.15 

0.05 

0.005 

0.001 

0.20 

• 0.08 

0.006 

0.0015 

0,30 

0.10 

0.008 

0.002 

0.40 

0.10 

0.01 

0.003 

0.50 

0.10 

0.015 

0.005 

0-.60 

0,15 

0.02 

0.00,8 

0.80 

0.20 

0.03 

0.01 

1.00 

0.10 



The optimization o£ the computational parameters was carried out in a 
separate study by considering 'the CO spectral lines in the region 2070- 
2220 cm“^ and by assuming a uniform CO distribution of one part per million 
by volume. 

Computations were made for four different values of the subinterval 
size- 6 . This ; is denoted in the computer program by' DEL (see Appendix C), 
Average transmlttances obtained for the entire range of CO fundamental 
band (between 2070-2220 cm”*^) are shown in Table 3.3 for the four subinterval, 
sizes. An attempt was also made to compute the transmittance for DEL =' 

3 cm~^ but was. eventually abandoned as the computational cost became top 
large. A value of DEL = 6 cm~^ was adapted for the present study simply 
because of the cost considerations. It was not considered necessary to go 
to even lower values of 6 in order to gain slight accuracy at higher cost. 

The number of adjacent subintervals (on each side of the main subinterval) 

\ 

contributing to the wing effects were also investigated. The results shown 
in Table 3.3 indicate that consideration of seven subintervals on each side 
yields sufficiently accurate results. It was not necessary to consider the 
influence of all neighboring subintervals in order to gain slight improvement 
(<0.02%) in the final result. 

Since averaging of the line widths over each decade in every subinterval 
is unavoidable (because of the approximations in the basic formulation of the 
band model) , no attempt was made -to consider the individual line widths even 
for CO and H 2 O lines. Average line widths for these molecules were found 
to have 10-15% standard deviation, ’The effect of using the average line 
widths on the overall transmittance is expected to be minimal. 



TABLE 3.3 


. Average CO Fundamental Band Transmittance as a 
Function of Various Quasi-Random Band Model Computational Parameters 


Parameters 

Range 

Average Transmittance 



(nondimensional) 

DEL (cm“ ) 

15 

0.81664 


10 

0.81523 


6 

0.81183 


5 

0.81050 

Number of 

4 

0.81225 

Intervals 

7 

0.81198 


8 

0.81194 


10 

0.81190 


all 

. 0.81183 



IV. MODEL CALCULATIONS 


In order to establish the feasibility of measuring the average concen- 
tration of pollutants • in the atmosphere (either from an aircraft mounted 
-instrument or from' -satellites) , it is essential to perform model calculations 
to evaluate upwellihg radiance and signal change under physically realistic 
conditions. By using the appropriate line parameters and employing the 
Lorentz' line-by-line -model and quasi-random band' model, upwelling radiance 
(and signal change) calculations were made in the spectral range of CO 
fundamental band (2070-2220 cm”^). Results were obtained to study the 
effects of different interfering molecules,, water vapor profiles, ground 
temperatures, and ground emittances on, the upwelling radiance and the signal 
change. In this section, the procedure for calculating the upwelling radiance 
and signal change is briefly discussed, information regarding data sources is 
provided, and finally,, the discussions of various results are presented. 

4.1 Procedure for Calculating the Upwelling Radiance and Signal Change 

In radiation modeling for pollution, measurement in a nonhomogeneous 
atmosphere, the upwelling radiation is calculated by dividing the atmosphere 
into an appropriate number of sublayers. Each sublayer is assumed to be 
homogeneous in species concentration, temperature and pressure. 

In a specified spectral interval in which a particular pollutant 
absorbs, the' total energy emergent from the atmosphere is obtained from 
Eq.(2,15). If in this interval, n independent measurements (corresponding 

to the number of homogeneous layers) could be made to find ^ 02 » ^Dn* 

then the uniform concentration of the pollutant in each layer (and therefore 
the concentration profile in the actual atmosphere) could be determined from 
Eq. (2,15). Because of low concentrations of pollutants in the atmosphere, 
however, n such measurements are not feasible. Thus, only one independent 
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measurement is usually made and an average value of the particular pollutant 
concentration in the atmosphere is obtained. Even if only one value of the 
pollutant concentration can be obtained from an independent measurement, it 
is essential to divide the nonhomogeneous atmosphere into several homogeneous 
layers for the purpose of data reduction. This is because the pressure, 
temperature, and amount of interfering molecules vary in the atmosphere, and 
spectroscopic parameters and pressure path lengths are strong functions of 
these variables# 

By employing the Lorentz line-by-line model for atmospheric transmittance, 
the upwelling radiance at the top of the atmosphere is obtained from Eq,(2.15) 
for each narrow spectral interval Aw . The exact procedure for doing 
this is to evaluate the average value of the Planch function for this interval 
first, then by using the mean value of the transmittance for the interval, 
Eq,(3.10), evaluate the upwelling radiance at the top of the atmosphere. The 
total upwelling radiance (E = Z at the top of the atmosphere for the 

entire spectral range A is obtained by summing the radiances of individual 
intervals. As pointed out earlier, for the present study the spectral range 
of interest is the range of CO fundamental band. The computer program for 
calculating the upwelling radiance in this spectral range is given in Appendix B, 
By employing the quasi-random band model, the total upwelling radiance 
at the top of the atmosphere can be evaluated from Eq#(2.15) in the following 
two ways , 

(a) First calculate the integrated Planck function for the entire spectral 
range A • Then by using the value of average transmittance for the 
range A , obtain the net (integrated) radiance at the top of the 
atmosphere. 

(b) Obtain the net (integrated) radiance at the top of the atmosphere for 
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each sub interval by calculating the Planck 'function and the average 
transmittance for that subinterval. Obtain -the total radiance for the 
entire range A by summing the integrated radiances of each sub interval. 

For the spectral range of CO fundamental band, total radiances at the 
top of the atmosphere were obtained by employing the above two procedures 
and these are 2,8102 E-05 and 2,8423 E-05 respectively. Procedure (b) was 
adapted for the present work because of its correctness in the formulation 
logic. Also, results obtained by this procedure were found to be in good 
agreement with the line-by-line results, A further simplification was 
introduced in procedure (b) in the final version of the computer program. 

Planck functions were calculated at the center frequency of each sub interval. 
This value was used as a mean for the sub interval rather than calculating it 
by integration over each 6 , No difference in the final results was observed. 
This is because Planck function is a slow varying function of the frequency 
(especially within the spectral range of subinterval 6 ). 

The signal change SC - AE (in watts/cm^-sr) can be calculated by 
employing Eq,(2.15) as 

SC = AE = I [E(w,t ) - E(w,t )]d03 , (4,1) 

•'Aw ° P 

where represents the transmittance of a "clean” atmosphere in which the 

pollutant concentration is zero, and refers to the transmittance of the 

atmosphere in the presence of the pollutant. The numerical procedure for 
evaluating Eq.(4.1) is identical to that described for calculating the 
upwelling atmospheric radiance. 

4.2 Data and Data Sources 

Upwelling radiances and signal change values were calculated in the 
spectral range of CO fundamental band. Contributions from the thermal 
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radiation and reflected sunlight terms in Eqs,(2,15) and (4.1) were 
taken into consideration while that from scattering terms were neglected. 

For the present study, the top of the atmosphere was considered to he 
10 km which is approximately the top of the troposphere. For numerical' 
calculations, this nonhomogeneous atmosphere was divided in. 10 layers of 
equal thickness. The sun zenith angle was taken to he zero (overhead sun, 

0 = o) for all calculations. 

The data- used in the present work were obtained from a number of well 
known sources and are considered to be highly reliable. 

The atmospheric temperature and pressure profiles were adapted from 
the U, S. Standard Atmosphere, 1962 [17]. Distributions of infrared 
active species, such as H 2 O, CO 2 , N 2 O, and O 3 were taken from McClatchey 
et al. [18]. Rotational and vibrational partition functions, needed to 
account for the temperature dependence of line strengths, were obtained 
from McClatchey et al. [2], Table 4.1 shows the range of variation of some, 
of these parameters. Values of, the parameters listed in the table refer to 
the conditions at the boundaries of each layer. Mean values of the parameters 
in each layer (corresponding to the values at layer centers) were used in 
actual calculations. Carbon dioxide and nitrous oxide are assumed to be 
uniformly mixed in the atmosphere. For these gases, average values for 

• - f 

the concentration used in the calculations are: CO 2 = 330 ppmV , and 

N 2 O = 0.28 ppmV. 
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TABLE 4.1 


Distribution of Pressure, Temperature, Water 
Vapor, and Ozone in the Earth’s Troposphere 


Altitude 

(fan) 

Temperature 

(“K) 

[Ref .17] 

Pressure 

(mb) 

[Ref .17] 

Water Vapor 
(p-pmV) 
[Ref. 18] 

Ozone Cone 
(ppmV) 
[Ref. 18] 

0 

288.15 

1013.25 . 

7.756 E+3 

2.663 E-2 

1 

281.65 

898.76 

6.068 E+3 

2.936 E-2 

2 

275.15 

795.01 

4.637 E+3 

3.239 E-2 

3 

268.66 

701.21 

3,187 E+3 

3.321 E-2 

4 

262.17 

616.60 

2.162 E+3 

3.391 E-2 

5 

255.68 

540.68 

1.399 E+3 

3.691 E-2 

6 

249.19 

472.18 

9.268 E+2 

4.118 E-2 

7 

242.70 

411.05 

5.731 E+2 

4.914 E-2 

8 

236.22 

356.52 

3,676 E+2 

5.973 E-2 

9 

229.73 

308.01 

1.586 E+2 

9.181 E-2 

10 

223.25 

265.00 

7.008 E+1 

1.315 E-1 


Science Applications Incorporated (NASA-Contractor , responsible for the 
development of the non-dispersive correlation instrument for pollution 
measurement [3,19]) has compiled spectral line parameters (position, strength, 
width, ahd lower energy level) for lines of CO fundamental band and for lines 
of other molecules which interfere with the CO band. In the calculation of 
transmittances, these line parameters are directly read from a tape provided 
by SAl [19] to NASA-Langley. 

Solar irradiances at the top of the atmosphere, at a few selected wave 
numbers, are also available from the SAI tape [19]. Values obtained from the 
tape, for the spectral range of present interest, are listed in Table 4.2. 

In the evaluation of contribution of the reflected solar radiation to the 
upwelling radiance, the solar irradiance for each spectral subinterval is 
obtained from a linear interpolation of values given in Table 4.2, 
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TABLE 4.2 

Solar Irradiances at the Top of the Atmosphere 
at Selected Wave Numbers 


Wave Number (cmT^) 


Solar Irradiance (watts cm“^ sr“^ (cm~^)“^) 


2061,86 

2105.26 

2150.54 

2197.80 

2247.19 


0.307 E-6 
0.323 E-6 
0.330 E-6 
0.343 E-6 
0.359 E-6 


4,3 Results of Model Calculations 

By employing the Lorentz line-by-line and quasi-random band model for 
atmospheric transmittance, upwelling radiance and signal change were calcu- 
lated for several illustrative cases. As indicated before, all calculations 
were made for the spectral range of CO fundamental band. 

The results of integrated upwelling radiance at the top of the tropo- 
sphere (i.e., at 10 km.) for different CO concentrations (uniformly distributed 
through the troposphere), in the presence of ^various interfering molecules, 
are illustrated in Fig. 4.1. The solid curves represent the results of the 
line-by-line model and broken curves for the quasi-random band model. As would 
be expected, the upwelling radiance E decreases with increasing CO concen- 
tration and with the inclusion of different interfering molecules. Inclusion 
of 0^ causes a slight decrease in radiance (not exceeding 0.5%) and it was 
difficult to illustrate this decrease in Fig. 4.1. The agreement between the 
line-by-line and the quasi-random band model results is seen to be excellent 
for the case of CO + H 2 O. The slightly lower radiances for the next two cases 
is attributed to the over-estimation of absorption by the band model. The 
reason for this lies in the assumption of random distribution of many lines 



(in the presence of interfering molecules) in the subintervals of the band' 
model. In the actual spectra, however, the lines are more closely spaced 
in some regions than in others. The reason for slightly higher radiance 
values, with the band model for the case of pure CO atmosphere (where 
relatively small number of lines are present), is not clear at this time. 

The variation in the signal change, AE , with the CO concentration is 
illustrated in- Fig. 4.2. These results follow the general trend of the 
results presented in Fig. 4.1# 

The influence of different amounts of water vapor on the upwelling 
radiance and the signal change is shown in Figs. 4.3 and 4.4 respectively. 
Increased water vapor concentration results in increased absorption in the 
atmosphere. This, in turn, results in lower values for upwelling radiance 
and signal change* It should be noted that the effect of CO concentration on 
the signal change would be relatively small in the presence of larger quantity 
of water vapor. 

Figures 4,5 and 4,6 show the upwelling radiances and signal changes for 
surface temperatures of 280, 290, and 300 ®K and a surface emittance of 0,8. 
The strong dependence of the upwelling radiance and signal change on the 
surface temperature is obvious from these results. The relatively lower 
radiance values obtained with the band model are indicative of slight over- 
estimation of absorption by this model. 

Figure 4.7 shows the variation of upwelling radiance for three different 
values of surface emittance and for a surface temperature of 288 As 

explained earlier, the radiances obtained from the band model are .lower than 
the line-by-line model because of overestimation of absorption by the band 
model. The relative increase of the difference for the lower e-values is due 
to lower total emission from the earth for the small values of surface 
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Fig, 4.2 Signal change as a function of CO concentration in the presence 
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Fig. 4.3 Upwelling radiance as a- function of CO concentration for different 
water vapor profiles, T = 288°K , e = 0.8 . 
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Fig, 4,4 Effects of water vapor concentration on the signal change 
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Fig. 4.6 Signal change as a function of CO concentration for 
.three different surface temperatures (surface 
emittance = 0.8). 
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emittance. In these 'cases, therefore, the increased absorption by the band 
model has a greater relative effect on the radiance* Figure 4.8 shows the 
signal change for the cases considered in Fig. 4.7 and the results indicate 
the same general trend. 

Figure 4.9 shows the variation of the upwelling radiance with the surface 
temperature for a fixed concentration of CO (1 ppm by volume) in the atmosphere 
and for e = 0.8. The strong dependence of radiance on the surface temperature 
may be easily explained on the basis of the Stefan^s law. However, because of 
the interference from the infrared active atmospheric molecules, the results 
obtained here do not exhibit an exact fourth power relationship. 

Figure 4.10 shows the variation of radiance with the surface emittance 
for a fixed CO concentration (1 ppm by volume) and = 288 As would be 

expected, the results indicate the linear dependence ‘of radiance on the surface 
emittance. Similar results were obtained in reference [3] for different values 
of the surface temperature and the sun zenith angle. In general, the ground 
emittance varies with the wave number. However, for the spectral range of 
CO fundamental band, it was shown in- reference [20] that the radiance is not 
influenced by a significant amount when the wave number dependent ground 
emittance is replaced by an averaged value. 

Figure 4.11 shows a comparison of the results obtained from the present 
line-by-line program (given the name LINBLIN) and another line-by-line program 
(called POLAYER) developed by the Science Applications Inc. [19] . In computing 
the total absorption coefficient at any wave number, the present program con- 
siders contributions from all the lines up to a fixed wave number location 
(on both sides) of 45.5 cm^^ from the wave number under consideration. This 
value of 45.5 cm”^ for the so-called wing effect was chosen after several 
numerical experimentations. The POLAYER, on the other hand, considers the 
effect of a fixed ntimber of 20 lines on each side of the wave number under 
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consideration* This causes the range of wing effect to change depending 
upon the density of lines in the spectrum* Thus, in some cases, POLAYER 
will not consider the influence of lines which are only 1 cm ^ away from 
the wave number under consideration* This, of course, will result in under- 
estimation of absorption* This, at least in part, is responsible for the 
higher integrated radiance obtained from the POLAYER program. The difference 
between the signal change curves , shown in Fig* 4.12, follows from the diff- 
erence in the radiance values. 
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Fig. 4.8 Signal change as a function of CO concentration for three' 
different surface emittances, T = 288“K . 
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Fig. 4.11 Comparison of upwelling radiance obtained from 
POLAYER program and the present LINBLIN program, 

T = 288 ®K , e = 0,8, On the average, LINBLIN 
values are 2% lower than POLAYER values. 
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Fig. 4.12 Comparison of signal change results obtained from the 
program POLAYER and the present LINBLlk program, 

T = 288°K, e = 0.8. On the average, LINBLIN values 
are 2% lower than POLAYER values. 



V. CONCLUSIONS 


A computer program has been developed for calculating the atmospheric 
transmittance and upwelling radiance in the spectral interval 2070-2220 cm“^ 
which employs the Lorentz line-by-line model. The spectral region selected 
in the region of CO fundamental band. The program allows for variation in 
CO concentration and accounts for the radiative contributions from the inter- 
fering gases such as H 2 O, CO 2 , N 2 O, and O^. The program has been optimized 
for accuracy and ininning time by varying such computation parameters as the 
interval width, sub interval width and the width of the region from which lines 
make significant contribution to the absorption at a certain frequency. 
Furthermore, this program has been diversified so that radiances can be 
computed for several values of surface temperature or surface emittance in a 
single run, thereby drastically reducing the computational cost. This line- 
by-line computer program is given the name LINBLIN, Upwelling radiance results 
obtained by using the present LINBLIN program are compared with the results ‘of 
SAI-POLAYER program [19]. The agreement between the two results is within 2%. 
For homogeneous gas columns, the results obtained by LINBLIN program compare 
well with available experimental results. This LINBLIN prograp can easily 
be extended to any line model and for any spectral region, 

A computer program similar to the line-by-line program has been developed 
xjhich employs the quasi-random narrow band model. This program was given the 
name QRANDOM, Comparison of results, obtained from LINBLIN and QRANDOM programs 
indicate small differences. 

By employing the LINBLIN and QRANDOM programs, model calculations were 
made to study the effects of different interfering molecules, water vapor 
profiles, ground temperatures, and ground emittances on the upwelling radiance 
and signal change. Physically realistic values of various parameters were 



used in the model calculations. This information is very useful in 
establishing the feasibility of measuring the concentration of pollutants 
in the atmosphere from a gas filter correlation instrument flown on an 
aircraft or mounted on a satellite. 
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APPENDIX A 

EXPLANATION OF SYMBOLS USED IN COMPUTER PROGRAMS 


A-i. Symbols Used in Program LINBLIN 


ALC ) 

ALV f 

ALX 

ALY 

ALl, AL2, 
AL3, AL4, > 
ALD ) 

DEL 

DLI-M 


EMI 


ELI, EL2, 
EL3, EL4, 
ELD 

FRL 

FRU 


FRB 


FRC 

FRl, FR2, 
FR3, FR4, 
FED 


GAM 


GR 


Individual line widths for CO and H 2 O lines 
respectively, cm~^ . 

Width of the narrow subinterval on each side 
of the line center , cm“ ^ . 

2 X (ALX) 


Average line width for the molecules H 2 O, C02» 
N 2 O, O 2 and CO respectively, cm“^ . 

Width of an interval, cm“^. 

Wave number region around an interval from which 
the contribution to the absorption coefficient is 
considered direct, cm~^. 

Surface emittance. 


Energies of the lower states for the lines of 
the molecules, cra^^. 

Lower frequency limit of the range (2070 cm"^. 
Upper frequency limit of the range (2220 cm“^. 
Wave number at the interval boundaries, cm**^. 
Wave number at the interval centers, cm“^. 

Wave numbers of the lines of the molecules, cmT^. 


1 + f(6) where £(6) = sec 0 for 
0 ^ 60° and f(0) = Ch (0) for 0 > 60° . 
Ch(0) is the Chapman function. 

Gradients used in calculation of HL 


A 

The same order subscripts 1,2, 3, 4, and D has been used for these 
molecules throughout the programs. 
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HS 

HL 


.INI, IN2, 
INS, IN4., 
IND 


Tabulated values of the solar irradiance at the 
top of the atmosphere, watts cm“^ sr~^ (cra“^)~^« 

Computed solar irradiance at the top of the ' 
atmosphere in an interval, watts cm“^ sr~^. 

An integer for the molecule which is 0 when 
average line width is used and 1 when individual 
line widths are used. 


JQ 


Number o.f different CO concentrations used in the 
computation ( 8 ) , 


KR 


Number of spectral intervals (150) * 


LC 


Number of layers into which the atmosphere is 
divided ( 10 ) . 


•LEI, LE2,. 
LE3,- LE4, 
LED 


Number of lines in the spectra of the 
molecules . 


NG 


Number of interfering gases. 


NGl, NG2, ) 
NG3, -NG4 / 

PCK 

PNTP 

PREC 

QV 

QVD 

RADNC 

RADNCE 


Identifying integers for the interfering molecules 
H 2 O, CO 2 , N 2 O and O 3 respectively, 

Planck’s function. 

Pressure at NTP, mbar 

Pressure at the center of the layers, mbar 

Altitude distribution of the molecules H 2 O, CO 2 , 
N 2 O and 0 ^ respectively. 

Concentration of CO, ppm by volume 

Radiance in each interval, watts cm“^ sr“^ , 

Integrated radiance, watts cm~^ sr“^. 


RCOM 


Radiance component due to reflected solar radiation. 


RPl, 

RP3, 

EPD 

SIl, 

SIS, 

SID 



Exponent to account for the temperature dependence 
of the rotational partition function for the molecules. 


Integrated intensities for the lines of the 
molecules , cm“^atom”^ . 
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TEMC 

TEMR 

TEtS 

TRC 

raTP 

TEA 

TRF 

TRM 

TRT 

TEX 

VPE 

WLIH 

WLl, 

it^N 

XLl, 


Temperatures at the layer centers, “K 

Reference temperature for line parameters, 

Surface temperature, “K 

Thickness of each layer (1 km) 

Temperature at NTP (273 ®K) 

Transmittances for different intervals and 
and altitudes for one gas. 

Combined transmittances of all the interfering 
gases. 

Average transmittance over the entire frequency 
range for each CO distribution. 

Transmittance for each interval and CO concentration, 
between the top of the atmosphere and surface. 

Transmittance for an interval and altitude including 
the contribution from CO. 

Vibrational partition functions for H 2 O, C 02 > N 2 O 
and O 3 . 

Wave number range of the x^ing contribution, cm“^. 

Weight factors for the 4-point Gauss-Leg endre 
quadrature formula. 

Wave numbers at which the solar irradiance at the 
top of the atmosphere is tabulated, 

XL2 Abscissa values for the 4-point Gauss-Legendre 

quadrature formula. 
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A"2. Symbols Used in Subroutine TRAILS 

AC Total absorption coefficient at the frequency 

specified by the subscripts, 

ACD Direct contribution to the total absorption 

coefficient, 

ACW Wing contribution to the total absorption 

coefficient, 

ACB Absorption coefficient at the lower boundary 

of the interval due to wing contributions, 

ACE Same as ACB, at the upper boundary of the interval, 

ACM Same as ACB, at the center of the interval, 

AL Average width for the lines of the molecule 

under consideration, 

ALA ^ Altitude dependent average width of the lines 

of a molecule^ 

ALB Altitude dependent individual widths of the lines of 

a molecule (for CO and H 2 O) . 

EL Energies of lower states of lines of the molecule 

under consideration. 

FACT Factor used in the computation of altitude dependence 

of integrated intensity, 

FR Frequencies of the lines of the molecule under 

consideration, 

FM Frequencies of the lines falling with an interval, 

FRG Frequencies of all the Gauss-Legendre quadrature 

points within an interval, 

FRS Frequencies at the subinterval boundaries within 

an interval, 

IG Number of Gauss-Leg endre points within an interval, 

IN An integer for each molecule which is zero if 

average line widths are used and 1 if individual 
widths are used for that molecule, 

LE Number of lines for the molecule under consideration. 
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I 


Number of sub intervals within an Interval* 

PART An altitude dependent factor which takes into 

account the vibrational and rotational partition 
functions* 

Pli Optical path length at the frequency under consid-- 

eration* 

f 

RP Exponent which accounts for the temperature 

dependence of the rotational partition function* 

SI Integrated intensities of the lines of the molecule 

under consideration* 

SXA Altitude dependent integrated intensities of lines 

of the molecule under consideration* 

SLl Gradient used in the calculation of the wing 

contribution between ACB and ACM* 

SL2 Same as SLl, between ACM and ACE. 

TR Transmittance at the frequency under consideration. 

WIDF Factor which takes into account the altitude 

dependence of the xc^idths of the lines * 
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A-3 . Symbols used .in PROGRAM QRANDOM 
TBA 


TRP 

TRIM 

TR6 

XI, T1 
X2, T2 


Combined transmittance of all the interfering 
gases. 

Combined transmittance including CO contribution 
for all CO distributions. 

Average’ transmittances for the entire frequency 
range for different altitudes and CO contributions. 

Transmittances for a single gas returned- by the 
subroutine. 

Convolution parameters for integration of direct 
contribution . 

Convolution parameters for the integration of wing 
contribution . 


As far as* possible, same symbols have been used in the programs 
LINBLIN and QRANDOM, Only the symbols which are defined differently 
for QRANDOM are explained here. 
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A-4. Symbols Used in SUBROUTINE SMITH 

AVSI Average value of intensity for the lines in one 

decade in an interval 

BIGI Intensity of the strongest line in an interval, 

cm.-^ atm“^, 

BIG Intensity values separating the five decades in 

each interval, 

ELE Altitude dependent energies of the lower states 

for lines in an interval. 

JB Number of adjacent intervals from which the wing 

contribution is considered. 

NSI Number of lines in a decade within an interval. 

SEE Intensities of the lines within a. decade in an 

interval. 

SIE Intensities of the lines within an interval, 

SLE Altitude dependent intensities of the lines 

within an interval. 

SSI Sum of the intensities of the lines in a decade 

within an interval. 

SUMSI Optical thickness between the top of the atmosphere 

and the altitude under consideration. 
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APPENDIX B 

LINE-BY-LINE COMPUTER PROGRAM TO- CALCULATE 
ATMOSPHERIC TRANSMITTANCE AND UEWELLING RADIANCE 

,B— 1 . Computer Program LINBLIN 


C 

C 


D I MENS I ON 
DIMENSION 
DIMENSION 

dimension 

DIMENSION 


ORIGINAL PAGE IS 
OP POOR QUALITY 


program L I NBL I N ( I NPUT t OUTPUT ) 

DIMENSION FRl nOO ) ,SI 1 ( 100 ) »£LX ( 100) 

FR2(4O0 } .512(400) .ELH(400 ) 

FR3 (750 ) .313(750 ) .EL3(750) 

FR4(1220) .S'14(1220).EL4(1220) 

FRD ( lOO ) .SI0( 100 ) •*ELO( 100) 

QVD ( 1 0 ) » QV'( 10.5). ALC ( TOO ) ♦ TRF (11.150). TRX (11.150). 

/TRT (150.8). PCK (11. 1'50 ) . GR ( 5 ) . RADNC ( 1 50 1 8 ) » TRM (10). RADNCE (10)'. 

/ VPF ( 1 0 . 4 ) » ttiN ( 5 ) . HS ( 5 ) . HL ( I 50 ) 

INTEGER G.X 

COMMON FRC (150). FRB (151) .PREC (10). TEMC ( 1 0 ) . QG ( 1 0 ) . TRA (11.150). 

/TEMR . XL i 1 XL2 . WL 1 » WL2-. ALX . LC . KR . PNTP . TMTP . DEL A . DL I M . WL IM . THC . 

/ALV( lOO )-.VP( lO) .ALY.DEL.PART( 10) .FACT(10).WlDF( lO> 

READ 10. FRL.FRU,.D£L»THC»PNTP»TNTP.TEMR.TEMS»EMI .ZEN 

READ 11. LC.JO.NG.LEl .LE2.LE3.LE4.LED 

read 10. DLIM.WLIM 

read 10’ '(OVD'( J) . J=1 .JQ) 

READ 15-» ( (QV(L.G).L-I .LC) »G=T'.NG) 

READ 10. {PREC(L ).L=1 »LC) 
read 10. (TEMC(L) .'L=4.LO 
read 12. XL1.XL2.WL1.WL2 
read 13. ALl .AL2 . AL3'.AL4 . ALD.ALX 
READ 10. RPl .RP2.RP3.RD4.RPD 
READ 17. ( (VPF (L.G) .L=l »LC) »G=1 .NG) 

■READ 11. NG1.NG2.NG3.NG4 
READ 11. INI ♦ 1N2. IN3. IN4, INO 
READ 18. (WN(N) .N^l .5) . ( HS (N ) . N= 1 . 5 ) 

READ 16. (FRl (X ) »SI 1 (X ) .ELI (X) i ALV(X > iX=l .LEI ) 

READ 14. (FR2(X ) .SI2 (X ) .EL2 (X) .X=l .LE2 ) 

READ 14. (FR3(X).Sf3(X).EL3(X).X=l.LE3) 

READ 14. (FR4 (X ) »SI 4 (X ) .EL4 (X ) .X=l .LE4 ) 

■ READ 16. (FRD(X) »SID-(X) .ELD(X).ALC(X).X = 1 .LED) 

10 FORMAT! iOF8. 2) 

11 FORMAT! 1615) 

12 FORMAT ( SF 10.6 ) 

13 format ( 10F8.3) 

14 FORMAT (2 (FI 1 . 3 .£ 1 3 . 4 .F 1 1 . 3 . SX ) ) 

15 FORMAT (8E1 0.3/2E1 0,3 ) 

16 FORMAT (2 (FI 1.3.E13.4.F11 .3..FS.3) ) 

1-7 FORMAT! 10f8. 4 ) 

IB F0RMAT(5F8.2 .SE8.3) 

CALCULATES THE NUMBER OF I NTER VALS-OEF I NES FREQUENCIES AT 

boundaries and centers 

• L8=LC+l 


THE .INTERVAL 
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DELA=0,5»DEL 

RK = ( FRU-FRL ) /DEL+0 * 1 

KR = RK 

, FR8( 2 )=FRL 
DO lOl K=1»KR 
FRB(K+1 )=FRB(K)+DEL 

101 FRC(K)=FR8(K)+DELA 
ALY=2**-ALX 

C INITIALIZES THE COMBINED TRANSMT TTAnCES OF ALL THE INTERFERING GASES 
DO 102 K=1»KR 
DO 102 L=1 *L8 

1 02 TRF ( L » K) = 1 . 

C DECIDE IF WATER VAPOR IS TO BE TAKEN INTO CONSIDERATION, IF SO» THE 
WATER VAPOR TRANSM 1 TTANCES RETURNED BY THE SUBROUTINE ARE COMBINED 
WITH TRF 

IF (NGl-2) 118,119,118 

119 CONTINUE 

• DO 1 07 L=1 ,LC 
VP(L)=VPF(L, 1 ) 

107 QG (L')=QV (L, U/1245* 

CALL T-RANS (FRl ,SI 1 ,ELl ,LE1 ,AL1 ,RP1 , INI ) 

DO 103 K=1 ,KR 
DO 1 03 L=1 ,LB 

103 TRF (L,!< )=TRF<L,K)*TRA(L,K-) 

C DECIDES IF CARBON DIOXIDE IS TO BE TAKEN INTO CONSIDERATION 
118 -IF {NG2-12) 165, 166,165 

166 CONTINUE 

DO 108 L=1 ,LC 
VPTL )=VPF (L,2) 

108 'QG(L')=QV(L,2) 

call TRANS ( FR2 , S I 2 , EL2 , LE2 , AL2 , RP2 , 1 N2 ) 

DO 104 K=1 ,KR 
DO 1 04 L=1 ,LB 

104 TRF(L,K)=TRF(L,K)*TRA(L,K) 

C DECIDES IF NITROUS OXIDE IS TO BE ‘ TAKEN INTO CONSIDERATION 
165 IF {NG3-7) 167,168,167 

168 CONTINUE- 

DO 1 0.9 L=1 ,LC 
VP (L )=VPF (L,3) 

109 QG<L)=Q.V(L,3) 

CALL TRANS(FR3,SI3,EL3,LE3, AL3,RP3, IN3 ) 

DO 105 K=1 ,KR 
DO 105 L=1,L8 

105 TRF(L,K)=TRF(L,K)»TRA(L,K> 


ORIGINAi; PAGE IS. 
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C DECIDES IF OZONE 15 TO' BE TAKEN INTO CONSIDERATION 
167 IF (NG4-19) 169»170»169 

170 CONTINUE 

DO 110 L=1'»L.C 
VP(L)=VPF(L*4) 

110 0G{L)=QV(L»4) 

CALL TRANS(FR4*SI4«EL4«LE4»AL4»RP4» IN4) 

DO 1 06 K=1 tKR 
DO 106 L=1 «L8 

1 06 TRF ( L ♦ K ) =TRF'< L VK ) *TRA ( L * K > 

169 CONTINUE . . • 

C calculates PLANCK FUNCTIONS FOR DIFFERENT ALTITUDES AND INTERVALS 
CONS = l 8**6«625-» l .E-Oy 
CNST=6. 625*0. 3/1 .38 
DO 115 K=1 iXR 
RNUM=DEL*C0NS*FRC {K)**3 
EEX=CNST*FRC (K ) 

PCK ( 1 . K ) =RNUsM/ ( EXP ( EEX/TEMS ) - 1 . ) 

DO ITS L=1 .LC 

115 PCK.(L+1 »K >=RNUM/ (,EXP (EEX/TEMC(L ) )-l . ) 

C CALCULATES SOLAR IRRADIANCE AT THE' TOP OF THE ATMOSPHERE IN EACH INTERVA 
DO 173 N=1 «4 

173' GR(N)={.HS(N+I )-H5(N) )/(WN(N+l )-WN(N) ) 

DO 174 K=1»KR 
N = 0 

177 N = N4-1 

IF (FRC(K)-WN(N) ) 174.175.175 

175 IF (FRC (K)-WN(N+1 ) ) 176.177.177 

176 HL(K)=HS(N3+GR(N)*(FRC(K)-VJN{N) ) 

HL<K )=HL(K)*D£L 

174 CONTINUE 

ZEN=ZEN/57. 29578 
GAM^l .+1 ./COS (ZEN ) 

DO 172 X=1.LED 
172 ALV(X )=ALC(X ) 

DO 178 L=1 .LC 
1-78 VP(L)=VPF(L. I ) 

C calculates TRANSM 1 TTANCES FOR all CONCENTRATIONS OF CO ONE BY ONE- 
C COMBINES THEM WITH TRF AND CALCULATES RADIANCE AT THE TOP TAKING 
C reflected SOLAR COMPONENT INTO ACCOUNT 
DO 111 J=1 ♦ JO 

DO 112 L=1 .LC 
112 QG(L )=QVD( Jr 

CALL TRANS (FRD. SID. eld. led. ALD.RPO. I ND ) 
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DO 113 K=1 .KR 
DO 114 L=1-.LB 

1 14' TRX<L»K)=TRF{LtK)*TRA(L«K) 

113 TRT(K« J)=TRX( 1 *K) 

DO 116 K=1»KR 

RCOM= ( 1 • -EM I ) * ( COS ( ZEN.) ) *HL ( K ) * < TRT ( K » J ) ) **GAM . 
RADNC(K»J )=EMI-»PCK{ I «K)*TRT (K* J)*l .E-07+RC0M 
DO' 1 16 L=2*LB 

COMP=PCK (LiK-)*(TRX(L »K )-TRX (L-1 »K ) )*1 *£-07 
1 1 6 RADNC ( K » J ) =RADNC ( K * J ) +COMP 
TRM1J)=0* 

RADNCEl J)=0» 

DO 117 K=1 «KR 
TRM( J)=TRM( J )+TRT{K» J)/KR 
1 1 7 RADNCE ( J ) =RADNCE { J > +RADNC ( K « J ) 

111 CONTINUE 

PRINT 60t { (TRT{K»J) .0=1 , JQ) ,K=1 «KR) 

PRINT 61. { (RADNC (K .J) .J=l .JO > .K=l .KR) 

PRINT 62. (TRM (J ) »J= 1 »JQ ),. (RADNCE (J ) »J=1 .JO > 

60 FORMAT! 1 HI/ (8F15«5‘) ) 

61 FORMAT ( 1H1/(8E1 5,4 ) ) 

62 format ( IHl///-// (8FK5.5////8E15.4 ) ) 

STOP 

END 
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SUBROUTINE TRANS ( FR , 5 1 » EL « UE » AL » RP ♦ I N ) 

DIMENSION PR (1220 ) *SI ( 1220 ) ♦ EL ( 1220 ) »SIA( 1220 ) «FRE (-40 ) »FRS(50) » 
/FRG<200> »ACD(200) « ACW( 200 ) « AC ( 200 * 1 0 ) » TR ( 200 ♦ 11 ) ♦ALB(IOO) 

integer X 

COMMON FRC ( 1 50 ) t FR3 ( 1 5 1 )•« PREC 1 1 0 ) ♦ I'EMC ( 1 U ) , OG ( 1 U ) t I RA (1 1 * I 50 ) ♦ 

/TEMR » XL 1 < XL2 » WL 1 ♦ WL2 , AlX » LC i KR * PNTP » TNTP * DELA » DL I M » WL 1 M» THC » 

/ALV( lOO ) f VP( IC ) *ALY.DEL»PART( 10) »FACT( 10) .V^DFC iQ) 

CST=SQRT (TEMR)ZPNTP 

C DEFINES VARIABLES WHICH ACCOUNT FOR THE VIBRATIONAL AND ROTATIONAL 
C PARTITION FUNCTIONS AND THE ALTITUDES VARIATION OF THE WIDTH OF LINES 
DO 120 L=1»LC 

PART(L)=VP{L)*( rEMRXTEMC(L) ) **RP 

FACT(L) = 1 .439*(TEMC<L)-TEMR)/t I'EMC ( L )*1 EMR ) 

120 WIDF<L)=CST*PREC(L)/SQRT(TEMC(L) ) 

x=o • • • 

C ONE INTERVAL CONSIDERED AT A TIME - FINDS ALL THE LINES FALLING WITHIN 
C IT - DETERMINES SUB INTERVAL BOUNDARIES WITHIN THE INTERVAL 
DO 123 K=1,KR 
M=0 
MP=M 

126 IF (X-LE) 124,J.25«125 

124 X=X+1 

IF (FR (X )-FR3(K > ) 126*127,127 

127 IF (FR(X )-FRB(K+l ) ) 12B, 125,125 

.128 M=M+1 

FRE(M)=FR{X) 

MP = M 

GO TO 126 

125 X=X-1 
N=l' 

FRS(N)3FR8 CK ) 

' IF <MP) 129,129,130 
130 CONTINUE 

DO 131 M=1,MP 
D I F=FRE ( M ) -FRS ( N ) 

IF (DIF-ALX) 132,132,133 

133 IF (DIF-ALY) 134,134,135 
1 32 FRS (N+ 1 > =FRE ( M ) 

N=N+1 
GO TO 136 

134 FRSCN+l )=FRE(M>-ALX 
FRS(N+2)=FRE(M) 

N=N+2 
GO TO 136 
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135 FRS(N+1 )=FRE;{M)-ALY 
FRS ( N+2 ) =FRE ( M ) -ALX 
FRS (N+3 )=FRE (M) 

N=.N.+3 

136 IF (M-MP) 137»138.»138 

137 DIF=FRE(M+1 )-FR:=>(N) 
IF(DIF-ALX) I39«139tl4u 

140 IFCDI'F-ALY) 141»14l»r42 
139 GO TO 131 

141 FRSlN+l )=FRE(M)+ALX 
N=M+l 

GO TO 131 

142 FRS(N+1 )=FRE (M)+ALX 
FRS ( N+2 ) =FRE ( M ) +ALY 
N = N+2 

131 CONTINUE 

138 DIF =FRB-(K+1' )-FRE (M)' 

IF (OIF-ALX) 143, 1 43 » 144 

144 IF (OIF-ALY) 14b., 145,146 

143 FRSCNFl )=FRB(K+1 ) 

NP=N,+ 1 

GO TO 147 

145 -FRS<N+1- )=FRE (M)+ALX 
FRS (N+2 )=FR3 (K+1 ) 

NP=N+2 

GO TO 147 

I 46 FRS ( N+ 1) =FRE ( M ) +ALX 
FRS (N+2 )=FRE (M l+ALY 
frS(N+3)=FRB(K+1 ) 


NP=N+3 
GO TO 147 

129 FRS(N+1 )=FR8 (K+1 ) 

NP=N+1 

147 NQ=NP-1 

C DETERMINES THE FREQUENCIES OF THE GAUSS-LEGENORE QUADRATURE POINTS 
C V'ilTHIN A SUB interval 
DO 148 N=1 ,NQ 
VAR = 0.5^^ (FRS(N+l )-FRS(N ) ) 

CON=0,5*(FRS (N+,1 )+FRS(N) ) 

I =4-a- (,N-1 ) + l 
FRG ( I ) =CON-VAR*XL I 
FRG( I+l ) =C0N-VAR*XL2 
FRG( 1+2 )=C0N+VAR*XL2 

148 FRG ( 1+3 )=C0N+VAR*XL1 
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IG=4*N.Q 
PI=3. 14159 

C EACH altitude CONSIDERED SEPARATEE Y-'o IRECT CONTRIBUTION EVALUATED 
C AT EACH GRID POINI AND W I NG CONlRIBUIION A1 THE BOUNDARIES AND 
C CENTER OF t Ht IN,. ERVAL 
DO 149 L=1«LC 
ALA=AL#WIDF(L) • 

ACB=0. 

ACM=0. 

ACE=0. 

DO 179 1 = 1 « IG 

ACD( I )=0« 

179 ACW( I )=0. 

IF (IN) 167»167»16S 
• 167 DO 150 X==l »LE 

DIF =A0S(FR(X)-FRC(K) > 

IF (OIF-WLIM) 151 ♦ 151 » 150 

151 IF (DIF-DLIM) 1 52 i 152 *153 * 

152 SI A(X ) =S I (X )«-PAR f iL) *EaPiEL ACT vL ) ) ' 

DO 154 1 = 1 * IG 

FD=FR (X )“FRG ( I ) 

DEN=P I (FD*FD+ALA*ALA ) 

154 ACD( I )=ACDi I ) +o I A v a ) *ALA/DEN 
GO TO 1 50 

153 S I A (X ) =S I (X ) *PAR f iL)->EaP i-EL k X ) A'C T ^ L ) } 

FB=FR (X )-FRB (K ) 

FM=FR (X }-FRC (K ) 

FE = FR (X )-FRB (K+1 ) 

PNUM=ST A { X ) *ALA 
ACB=ACB+PNUM/(p‘l'J‘Fd**2 ) 

ACM=ACM+PNUM/ (P I *FM*42 ) 

ACE=ACE+PNUM/ (P I *FE**2 ) 

150 CONTINUE 
GO TO 1 74 
168 DO 169 X=1 *LE 
, DlF=A8btFRlX)-FRC(K) ). 

IF (DIF-V/LIM) 17u.l7u*169 

170 IF (DIF-DLIN) 171*171*172 

171 SIA(X)=S1 (X)#PART{L)*EXP{EL(X)*FACT(L> ) 

ALB(X)=ALV(X)*WIDF(L) 

DO 173 1 = 1 * IG 
FD=FR(X )-FRG( I ) 

DEN=PI*- (FD*FD+ALB IX)**2) 

173 ACD( I )=ACD( I ) + S I A < X ) *ALB C X ) /DEN 
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GO TO 169 

172 SIA(X)=SI (X )*PART (L) *EXP (EL <X )*FAcr IL) ) 

ALB (X ) =ALy (X)*Wl'DF(L ) 

FB = FR (X )-FRB(K ) 

FM=FR(X )-FRC(K ) 

FE=fR C X ) -F-R8.(:K,+ 1 ) 

PNUM =S I A ( X ) *ALB ( X ) 

ACB= ACB+PNUM/ ( P I *FB**2 ) 

ACM=ACM+PNUM/ (.P [ «-FM**2 ) 

ACE=ACE+PNUM/ (PT *FE**H ) 

169 CONTINUE 

C WING CONTRIBUTION EVALUATED AT EACH GR I-D POINT BY INTERPOLATION 

174 SL1 = (.ACM-ACB)/(FRC(K)-FRB(K) ) . - 

SL2= ( ACE-ACM )/ (FRStK+l )-FRC (K) ) 

DO 156 I=i ♦ IG 
DIF=FRG( I )-FRB.(K) 

IF (OIF-DELA) 157»158.158 

157 ACW( I )=ACB+SL1*DIF 

GO TO 156 • ' 

158 ACW( I )=ACM+bL2*(DIF-DELA) 

156 CONI INUE 

C total A650RPII0N C'OEFF I C I LN I'U bVALv^ATtU A, EACH GRIDPOINI 
DO -1'61 I=i».-I,G 

161 ACd vL)=ACD( n+ACW( I > 

149 CONTINUE 
LB=LC+l 

CONST=O.I*TNIP*iHC/PNlP 

C TRANSMITTANCE EVALUATED -AT EACH GRID POINT 
DO 163 I=1»IG 
PL=0. 

TR'( I .LB )=EXP (-PL ) 

DO 163 M=2.LB 
L=L8+1-M 

SPTR=PREC ( L ) *OG ( L ) *AC ( 1 . L ) / JEMC ( L ) 

PL=PL+CONST*SPTR 

163 TR ( I .L ) =EXP(-PL) 

-C AVERAGE TRANSM I l‘ I ANCE EvALuA I ED FOR EACH INTERVAL At’ EACH AL'i I i uDE . 
DO 164 L=1,.LB 
TRA(L»K)=0* 

DO 164 N=-l»NQ 
VAR=0*5*(FRS (N+1 )-FRS (N) ) 

I =4* (N-1 )+l 

• • SUM1=TR(I ,L)+IRC I+3.L) 

SUM2=TR (I+l.L)+lR(I+2.L) 

SUM=WL1 *SUM1 +WL2*SUM2 

1 64 TRA ( L » K ) =TRA ( L . K ) +SUM* VAR/DEL 
123 CONTINUE 

. RETURN 
END 
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APPMDIX C 


QUASI-RANDOM COMPUTER PROGRAM TO CALCULATE 
ATMOSPHERIC TRANSMITTANCE AND UPWELLING RADIANCE 


Computer Program QRANDOM 


PROGRAM QRANDOM (INPUTS OUTPUT) 

INTEGER X»G»W 

O I MENS.I ON FR 1 ( l UO ) ♦ i, I 1 ( 1 uu ) . hL 1 l 1 uu ) 

DIMENSION FRaCAUO ) ,S 12(400 ) <ELlit4UO )' 

D I mens I ON FRD ( 1 00 ) . S I D ( 1 00 ) » ELD (100) 

■DIMENSION OVD( 10) t OV Uui 4 } » VPF ( lu«4 ) « wNt 5) i5)»GR(5) * . 

/TRA (.26 » 11 ) ♦ TRF { 5 1 » 1 1 ♦ 8 ) t HL ( 25 ) « PCK ( 25 . 1 1 ) » TRFM ( 1 1 ♦ 8 ) » TRT ( 25 » 8 ) t - 
/RADNC(25»8) »TRM{a,}tRADNCE(8) 

COMMON DEL* XI ( 26 ) » II (26) *X2 (21 ) * 12(21 ) *TRG(26* 1 1 ) *FRB127) *FRC.(26) » 
/PREC(IO) *TEMC(lO) . QG ( 1 0 ) * VP ( 1 0 ) ♦ PNTP ♦ J N I P t TEMR*LC *KR * JB , THC 
READ 10* FRL*FRU»DEL*THC,PNTP* INTP*TEMR*TEMS*EMI *ZEN 
read 11* LC* JQ»NG* JB *LE1 *LE2*LE3*LE4.LED 
READ lO* (QVD ( J ) ♦ J=1 * JQ) 
read 10. (PREC (L).*L=1 *LC) 
read 10* ( IEM'C(L) *L=1 *LC) 

read 12. ( IQV IL *G ) *L=1 *LC ) *G=1 *NG) 

READ 14-* ALl iAL 2,AL3*AL4*ALD 
READ 10* RPl .RP2.RP3 ,RP4-.RPD 

read 11* MP-1 *NR1 .MP2.NR2.MP3.NR3*MP4*NR4iMP0»NRD 
READ 1 5 * ( X 1 ( W ) * w = l * 26 ) 

READ 15* (Tl (W) *W=1 *26) 
read 15* (X2 (W) *.W=l *21 ) 

READ -15* (T2(V/) *W=l .21 ) 

read 15* ( (VPF(L*G) *L=1 *LC) *'G=1 *NG) 

read 11* NGl *NG2*NG3.NG4 

READ 18. (WN(N) *N = 1 *5) * (Hia(-N)*N=l .5) 

READ 17* (FRl (X ) *31 1 tX.) *EL l'(x ) »A=1 *LE1 ) 
read 17* (FR2(X ) *SI2 (X ) *EL2 (X) *X=1 *LE2: 

READ 17* (FRD(X ) .SIDCXO *ELD(X) *X=1 .LED) 

1 0 FORMA ( ( 1UF8.2) 
n FORMA I ( 1615 ) 

12 F0RMAT(8E10.3/2EIU*3) 

14 FORMAT (10F8 .3) 

15 FORMAT! 10F8. 4 ) 

17 FORMAT (2(F11«3.E13*4,F11.3.5X)) 

18 FORMAT (5F8 *2* 5E8»3) 

C CALCULAIEb I HE FREQuENCIEo At i HE IN,'ER./AL BOv^NDARIEo AND CENIERJ FOR 

C UNSHIFTED MESH AND 1N1IIALI,4E3 I HE REbUL I AN i l RANcsM I M Aj^CES 

RK= (FRU-FRL )/DEL+0, 1 
KR=RK 

0ELA = 0.5-i^DEL 
FR0 ( 1 ) =FRL 
DO 100 K=1 *KR 
FRB(K+1 )=FRB(K)+DEL 
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100 FRC(K)=FRB(K )+PELA 
LB=l_C+l 

DO 1-01 L=1»L8 
DO lOl K=ltKR 

1 0 1 TRA ( K » L ) = 1 . 

C OEC-l.DES WHICH GASES ARE TO BE TAKEN INTO CONSIDERATION IN THE ORDER 
C WATER VAPOR* CARBON DIOXIDE* NI fROUi OaIDE AND 0-^.ONE AND COMBINES 
C THEIR transmit I ANCEi> OBlAlNED FROM I HE Si^BROOTlNh wIIH , HE REoOLTANI 
C TRANSMITTANCE 

IF <NG1-2) !02*-103* 102. 

103 CONTINUE 

DO 104 L=1*EC 
QG(E)=QV<L* 1)/1245* 

104 VP<L)=VPF(L* 1 ) 

CALL SMITHIFRl *SI 1 , ELl *LE1 * ALl *RPl »'MP1 *NR1 ) 

DO 105 L=1*LB 
DO 1 05 K~1 ,KR 

105 TRA(K*L )=TRA (K»L)*TRG<K*L) 

102 IF (NG2-12') 106*107, 106 

107 CONJ.INUE 

DO lOa L=1 ,LC 
QG(L,)=QV(L*2 > 
lOB VP<L)=VPF{L*2) 

call SM I T H ( FR2 * S 1 2 , EL2 * LE2 * Al 2 * RP2 , MP2 . NR2 1 
00 l09 L=1 *L3 
DO 109 K=1*KR 

109 TRA (K,L )=TRA (K ,L) »-| -RG(K*L)- 

106 IF (NG3-7). 110,111,110 

111 CONTINUE 

DO 112 L=1*LC 
QG(L)=0V(L*3) 

112 VP(L >=VPF (L,3) 

CALL SMITH(FR3*SI3*EL3,LE3*AL3*RP3,MP3,NR3) 

DO 113 L=1 *L8 
DO 113 K-1 *KR 

113 TRA (K*L )=TRA {K,L)-»TRG(.<*L) 

110 IF {NG4-19) 114*115*114 

115 CONTINUE 

DO 116 L=1 *lC 
QG(L)=QV(L*4 ) 

116 VP (L )=VPF (L*4 ) 

CALL SM1THIFR4 *S I4*EL4*LE4* AL4*RP4*MP4*NR4) 

DO 117 L=1 *LB 
DO 117 K=1 ,KR 
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117 TRA(K»U)=TRA (K,l_)*TRG(K»L) 

114 CONTINUE 

DO 154 L=1»U8 
DO 154 K=1-,KR 
KK=2*K 

154 TRF(KK«I_. 1 )=TRA (K »L) 

DO 1 18 J=2t JO 
DO 119 L=1 «LC 
QG (L)=QVD( J) 

1 19 VP(L)=VPF(L» 1 ) 

call SMITH! FRO ♦ S I D t ELD » LEO » ALD» RPD * MPD ♦ NRO ) 

DO 120 L=1 »LB 
DO 120 K=I ,KR 
KK=2*K 

1 20 TRF ( KK » L t J ) =TRA ( K « L ) RG ( K ♦ L ) 

118 CONTINUE 

C calculates THE FREQUENCIES A I I HE IN’IERvAL BOUNDARIES AND CEN I ERb 
C FOR SHIFTED MEbH AND INDIALlZhS RtowL'l AN f VRANoM 1 T'l-ANLtS 
KR=KR+l 

FR6( 1 3=FRL-DELA 
DO 121 K=ltKR 
FR0(K+i )=FRB(K)+DEL 
121 FRC(K )=FRB IK >-rDELA 
DO 122 L=1 «LB 
DO 122 K=1 »KR 
. 122 TRA(K^L)=1. 

C CALCULAIEb I RANSM I I i ANCES AND COMBlNEo i HEM wlTH t Ht REoULfANi THE 
C SAME way as for uNi.HIFiED MFSH 
I- (NGl-2) 1 23 » 124 *123 

124 CONTINUE 

DO 125 L=1»LC 
QG(L)=QV(L« 1‘)/1245. 

125 VP (L )=VPF (Lt 1 ) 

call SM.ITH{'FR1 »SI 1 »EL1 »LE1 «AL1 «RP1 »MP1 iNRl ) 

DO 126 L=1 tLS 
DO 126 K= 1 .KR 

126 TRA(K»L) = TRA<.<»L)*lRGlK»L> 

123 IF (NG2-12) 127»128»127 

128 CONTINUE 

DO 129 L=1»LC 
OG(L )=QV(Li2 ) 

129 VP(L)=VPF(L»2) 

call SMITH! FR2 ♦ S I 2 . EL2 » LE2 » Al2 » RP2 . MP2 » NR2 ) 

DO 130 L=1«LB 
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DO 130 K=1 »KR 

130 TRA (K»L') = 1 RA lK»L)-«-|RGvK»L) 

127 IF <NG3-7) 131 » 132 ♦131 

132 CONTINUE 

DO 133 U=itLC 
QG(L)=QV(L*3) 

133 VP{L)=VPF{,Lt3) 

CALL SMITH (FR3»i>l3^EL3^LE3* AL3tRP3»MP3fNR3) 

DO 134 L=1 tLB 
DO 134 K=ltKR 

134 TRA(KtL)=TRA(K,L)*lRGtK»L) 

131 IF (NG4-19) 135»136^13b 

136 CONTINUE 

DO 137 L=ltLC 
QGCL.) =QV(L»4) 

137 VP (L ) = VPF (Lt4 ) 

CALL SMI I H(FR4,^SI4^EL4iLE4%AL4^RP4^MP4^NR4) 

DO 138 L=l^LB 
DO 138 K=1 fKR 

133 TRA(K^L)=TRAlKtL)*lRGtK^L) 

135 CONTINUE 

DO 1 55' L=1 ♦lb 
DO 1 55 K=1 ♦KR 
KK=2*K-1 

155 TRFIKK^L^ 1 )= IRA (K^L) 

DO 139 J=2^JQ 
DO 140 L=1^LC 
OG (L )=QVD ( J ) 

140 VP (L >=VPF (L^ 1 ) 

CALL SM I I H ( FRD ♦ I D ♦ ELD ♦ LED ♦ ALD^ RPU ♦ MPU ♦ NRD ) 

DO 141 L=1»LB 
'DO 14 1 .K=l *KR 
KK=2»K-1 

141 TRFIKK^L^ J ) = TRA IK^L)-»lRG(K^L) 

139 CONTINUE 

KR=KR-1 

C CALCULATES MEAN TRANSM I TTANCES BY AVERAGING RESULTS FROM SHIFTED AND 
C UNSHIFTED MEiH CALL»jLA I 1 ON J ANu FiNmLLV OoTAlNiNG i nc. AVc.RAGc: .'RSN_>M i T TAN 

C FOR the ENIIRE FREQutNC i KANGG 
DO 143 J=1 ♦JO 
DO 144 L=l ♦LB 
DO 145 K=i^KR 
KX=2*K-1 
KY=2*K 
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KZ=2*K+1 

1 45 TRf ( K ♦ L » J ) =.( TRF ( K X ♦ L * J ) +TRF ( K Y ♦ L » J ) +TRF ( KZ ♦ L * J ) ) /3 • 

TRFM(L»J')=0* 

DO 144 K=ltKR 

1 44 TRFM ( L » J ) = I RFM (L»J)+lRF(K»L*J) /'KR 
DO 143 K=1»KR 
143 TRT(K» J)=TRF(K* 1 »J) 

C CALCULAIES SOLAR IRRADIANCE A I' iHE I OP OF (HE A i,MO ^PHtRE FROM i HE 
C TABULAtED VALuto -tJ.T Ll NEAR' IN (tRPOLAT 1 ON iN tAEn INIeRVmL 
DO 1-46 N=U4 . 

146 GR(N)=(HS(N+1 )-HS.{N) )/(WN<N+l )-WN(N> ) 

DO 147 K=1 »KR 

N=rO 

150 N=N+1 

IF (FRCIK)-WN (N) ) 147»148»148 

148 IF (FRC IK )-WN (N+1 n 149,150*18^ 

149 HL(K) = OEL^t-<HSlN)+GRlN)*vFRC(K)-wN\N ) ) ) 

147 CONTINUE 

C CALCULATES PLANCK FUNCTION FOR EACH ALTITUDE AND FREQUENCY INTERVAL 
C0NS=18.-»-6.625*-l .E-U7 
CNST=6.625*u.3/l • JO 

DO 142 K=1*KR 
RNUM=DEL*CONS-)tFRC (K ) 

EEX=CNST*FRC(K ) 

PCK ( K * 1 ) =RNUM/ ( EXP I EEa/ I EM o ) - 1 . ) 

DO 14'2 L== l * LC 

PCK(K*L+1 )=RNUM/(EXP(EEX/IEMCIL) )-l • ) 

142 CONTINUE - ‘ . 

C CALCULAIES radiance Af IHE TOP OF i HE ATMOSPheRe ''INCL^^UING 1 HE 
C REFLECTED SOLAR COMPONENT 
ZEN=COS (ZEN/57.29578) 

GAM=1 .+1 ./ZEN 
DO 151 J=1»JQ 
DO 152 K=1,KR 

RCOM=( 1 .-EMI >*2EN*HL(K;)*(TRT(K, J)ji(-*GAM 
RADNC (K, J )=EMI*PCK (K , 1 )»TRT(k* J)*l .E-07+RC0M 
DO 1.52 L=2,LB 

COMP=PCK (K*L )■» (TRF(K,L*J)-TRF(K,L— 1 * J) )*1 .E-07 

152 RADNC(K» J.)=RADNC (K* J)+COMP 
TRM(J)=0. 

RAONCEt J)=0* 

DO 153 K=1 *KR 
TRM( J )=TRM( J )+TRT(K« J)/KR' 

1 53 RADNCE ( J ) =RADNCE ( J ) +RADNC ( K » J ) 

151 CONTINUE 

PRINT 60* ( (TRT(K.J) *J= 1 * JQ) *K=1 *KR> 

PRINT 62, ( CRADNC (K,J) *J=1 * JO )-«K = l ,KR). 

PRINT 60, ( (TRFM(L,J) . J=1 , JQ) *L=1 ,LB) 

■ PRINT 61* (TRMC J)*,J=1 ♦ JO) , (RADNCE( J) , J=1 ♦ JO) 

60 FORMAT! 1H1//(8F15. 5,/) ) 

61 F0RMAT( 1H1/////<SF15.5////8E15.4 ) ) 

62 FORMAT! 1H1//(8E15. 4,/) ) 

STOP 
END 
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SUBROUTINE SMITH(FR«S1 ♦ EL ♦ LE » AL » RP t.MP ♦ NR ) 

INTEGER XtW 

D I MENS I ON FR-( 1 220 ) » S I { 1 22-0 )•-. EL { 1 220 ) ♦ S I E C 200 ♦.26 ) ♦ SLE ( 200 ♦ 26 ) . 
/ELE(200^26) .SFE ( iOO ♦5i26) ♦8'IG'r(26 ) ♦ B IG ( 6 ^26 ) ♦ NS I ( 5 ♦26 ) ♦SSI (5^26) ♦ 

/ AVS-I { 5 ♦ 26 ) ♦ SUMS 1 ( 5 ♦ 26 ) ♦ PART ( 1 O ) ♦ FACT (1 O ) ♦ W t OF { 1 0 ) 

COMMON DEL^Xl (26) ♦Tl (26) iX2(21 )^T2(21 ) ♦TRG(26^ 11 ) ♦FR8(27) ♦FRC(26)t 
/PREC( 1.0 )-iTEMC( lO) ,QG( 10 ) >VP( 1 0 ) ♦ PNTP ♦TNTP , TEMR .LC ♦ KR ♦ JB ♦ THC ‘ 

C SETS INTEGRATED INTENSITIES AND ENERGIES FOR EACH LINE TO ZERO 
DO 120 K=1 ♦KR 
DO 120 ,M = 1 1 MP 
SIE(MtK)=0. 

120 ELE(M^K)=0. 

X=0 

C DISTRIBUTES THE LINES INTO FREQUENCY INTERVALS 
DO 121 K=1^KR 
M = 0 

124 -LF (X-LE) 1 22^1 23 ♦123 

122 X = X-tl 

IF (FR (X )"FRB (K ) ) 124^125^125 

125 IF (FR (X )-FR8 (K+1 ) ) 126^127^127 

126 M=M41 
SIETM^K )=S I CX) 

ELE(M^K)=EL(X) 

. GO TO i 24 

127 IF (K-KR) 128^r23^123 

128 X=X-1 

121 CONTINUE 

123 CONTINUE 

•C SETS OPT-ICAL DEPTH IN EACH INTERVAL AT THE TOP OF THE ATMOSPHERE TO ZER< 
DO 151 K=1,KR 
DO 151 I=l^5 

151 SUMSI ( I ♦K)=0* 

LB=LC+l 

PI=3..141592654 

DELA=0.5*DEL 

CONST=0 . 1 *TNTP*THC/PNTP 

CST=SQRT (TEMRl/PNTP 

ULB=0. 

BSU=0. 

C EACH layer is CONSIDERED SEPARATELY FROM TOP DOWN 
DO 150 MM=l^LC 
L=LC+1-MM 

PART(L )=VP(L.)*(TEMR/TEMC(L) ) **RP 

FACT (L ) = 1 »4 39* t TEMC (L )-TEMR )/ ( TEMC (L )*TEMR) 
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, W I DF { L > =CST*PREC ( L ) /SORT ( TEMC (L ) > 

•SPTR=PREC (L-)*QG(L)/TEMC(L) 

ALA = AL*vnOF(L) 

ULA=CONST*SPTR/ <P I*ALA,). 

ULB=Ul-8+ULA 

BSU=BS'J+ULA*ALA 

ALB=BSU/ULB 

RH0=ALB/DELA 

C CALCULATES THE ALTITUDE DEPENDENT INTENSITY 
DO 137 K=1.KR 
DO 137 M=1»MP 

137 SLE(M»!<)=SIE.(M,K)^PART(L)»EXP(’PA'CT(L)*ELE<M»K) ) 

DO 129 K=1«KR 
DO 129 1=1.5 

DO 129 N=1.NR 

129 SFE(N. I .!<)=0, 

C ASSIGNS LINES IN THE INTERVAL TO THE FIVE INTENSITY DECADES 
DO 130 K=l.KR 
6 IGI (K ) =SLE ( 1 .K ) 

DO 131 M=2»MP 

IF (B I G I (K ) -SLE (M .K ) ) 1 32 . 1 31 » 1 3-1 

132 B I G I ( K ) =SLE( M . K ) 

131 CONTINUE 

DO 133 1=1.6 
IX=-I+1 

133 BIG! I )=BIGI (K)*10.**IX 

DO 130 1=1.5 

N = 0 

DO 130 M=1.MP 

IF (SLE(M.K)-BIG( I .K) > 134.134.130' 

134 IF (SLE<M.K)“BIG( I+l .K) ) 130.130.135 

135 N=N+1 

SFE { N . I . K ) =SLE ( M , K ) 

130 CONTINUE 

-C CALCULATES AVERAGE I NTENS 1 TY' FOR .EACH DECADE 
DO 138 K=1.KR . .. 

DO 138 1=1.5 

NSr < I iK)=0 
SSI ( I .K)=0. 

DO 139 N=1 «NR 

IF (SFE(N.I.K)) 140.140.141 
141 NSI C I .K )=NSI ( I .K)+l 

139 SS I ( I . K ) =SS I { I . K ) +SFE < N . 1 ♦ K ) 

140 IF (NSld.K)) 142.142.136 
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142 NSI < I »K)=1 

136 AVSI ( I 1 K-)=SSI { I *k )/NSI U ,K) 

SUMS I C 1 «K)=SUMSi { I ♦Kl+AVSI t I *K)^«-ULA 
138 CONTINUE 

C calculates TRANSMITTANCE IN EACH, INTERVAL AT ALL ALTITUDES CONSIDERING 
C DIRECT AND WING CONTRIBUTION 
DO .143 K=l »KR 
FRUC=FRC(K) 

TRG(KiL)=l » 

DO 143 J=1VKR 
JA=I ABS ( J-K ) 

IF (JA-JB) 152* I 52 ♦153 
153 TRD=I . 

GO TO 143 
152 FRBI=FRB{J) 

- 2I=FRUC-FRBI“DELA 
EPSI=ZI/DELA 
TRD=1 . 

DO 144 I=I^5 

NSJ=NSI { I ♦ J ) 

XI=SUMSI{ I ♦ J) 

RES=0 . 

IF (J-K) 145»146fl45 

146 DO 147 W=1 426 

Y=EXP(-RHO*RHO#Xi/(Xl (W)»X1 ( ) +RHO*RHO ) ) 

147 RES=RES+Y*T1 ( W > 

GO TO 144 

145 DO 148 W=1 ♦Bl 

Y=EXP ( -RH0*RH0*X 1 / { EPS I -X2 ( W) 1 **2 ) 

148 RES = R£S-J-Y*T2(W> 

RdS=RES/6. 

144 TRD=TRD*RES4«-NSJ 

143 TRG(K»L }=TRG(KtL)*TRD 
150 CONTINUE 

C DEFINES TRANSMITTANCES AT THE TOP OF THE ATMOSPHERE AS UNITY 
DO 149 K=1 fKR 

149 TRG(KiLQ)=1. 

RETURN 

END 
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